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The propagation of light in a scattering medium is described as the motion of a special kind 
. of a Brownian particle on which the fluctuating forces act only perpendicular to its velocity. This 

enforces strictly and dynamically the constraint of constant speed of the photon in the medium. 
A Fokker-Planck equation is derived for the probability distribution in the phase space assuming 
the transverse fluctuating force to be a white noise. Analytic expressions for the moments of the 
displacement (x n ) along with an approximate expression for the marginal probability distribution 
function P(x,t) are obtained. Exact numerical solutions for the phase space probability distribution 
for various geometries are presented. The results show that the velocity distribution randomizes in 
a time of about eight times the mean free time (8t*) only after which the diffusion approximation 
becomes valid. This factor of eight is a well known experimental fact. A persistence exponent of 
0.435 ±0.005 is calculated for this process in two dimensions by studying the survival probability of 
the particle in a semi-infinite medium. The case of a stochastic amplifying medium is also discussed. 
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I. INTRODUCTION 



The propagation of light through a stochastic medium is traditionally described in the context of astrophysics by a 
Boltzmann transport equation for the specific intensity I(r, f2, t) in a heuristic Radiative transfer theory [jlj. However, 
since the general analytic solutions are unknown, one resorts to the diffusion approximation which can be shown to 
arise out of the Radiative transport equation in the limit of large length scales L >> I*, where I* is the transport mean 
'— '! free path of light in the medium [jyj2]. Recently there has been considerable interest in the description of multiple 
light scattering at small length scales (L ~ I*) and small time scales (t ~ t* where t* is the transport mean free 
time) , both from the point of fundamental physics (3) and from the point of medical imaging, where the early arriving 
'snake' photons are used to image through human tissues It has been experimentally shown that the diffusion 

■ approximation fails to describe phenomena at distances of L < 81* Moreover the diffusion approximation which is 
strictly a Wiener process for the spatial co-ordinates of a particle is physically unrealistic. It holds in the limit of the 
mean free path I* — > and the speed of propagation c — > 00 while keeping the diffusion coefficient D Q — ^- constant. 
Thus the diffusion approximation neither accounts for a finite mean free path nor for a finite and constant speed 
of the particle which is characteristic of light propagation in a stochastic medium. While approximately describing 
light as a particle the constancy of speed should be preserved at the very least. Hence it is of importance to develop 
better and alternative schemes to the diffusion approximation and also address the difficult question of the process of 
g ' randomization of a directional beam in such media. 

*^ ■ For a particle moving with fixed speed in a one- dimensional disordered medium, it has be shown that the probability 
distribution function P(x, t) for the displacement satisfies the telegrapher equation exactly. However, generalizations 
of the telegrapher equation to higher dimensions have been shown not to yield better results than the diffusion 
approximation Q. Recently there have been a few attempts to overcome the shortcomings of the diffusion approx- 
imations and attack this problem using the concept of photon paths. In Ref. ||, a Monte-Carlo approach was used 



to simulate photon paths and calculate their probabilities. An important advance was made in Refs. [|10 11 1, where 
the propagator for photons in highly forward scattering media was expressed as a Feynman path integral. However, 
this attempt has had only limited success in that it was possible to calculate the probability distribution subject to 
the constraint of constant photon speed only in the weaker (average) sense i.e., L [(^-) 2 — c 2 ] dt = 0. Moreover in 
addressing the backscattering from a semi-infinite medium |l0[ and reflection/transmission from a finite slab JTl| , the 
absorbing boundary conditions have not been rigorously implemented and it would be inappropriate to compare these 
to experimental data. It should be mentioned here that the Ornstein-Uhlenbeck (O-U) theory of Brownian motion 
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|L2| would also be able to incorporate the finiteness of the mean free path and a well defined root-mean-squared 
(rms) velocity but assuming, of course a distribution of speeds. This process has been compared with Monte-Carlo 
simulations and used to explain the lowering of the effective diffusion coefficient measured in pulse transmission 
experiments through thin slabs [p|. It can be shown that the finite r.m.s speed defined by the fluctuation-dissipation 
theorem for the O-U process is a stronger global constraint than the average speed constraint imposed in Ref . JlO|,|ll| . 

The next important step in describing these photon random walks with a constant speed was undertaken in Ref. 
fL4j , where the authors describe this process as a non- Euclidean diffusion on the velocity sphere and intuitively put 
down a kind of a general Boltzmann equation for photons in a highly forward scattering medium. The solution to 
this equation was expressed as a path integral, which was then evaluated by a standard cumulant decomposition |l5| l 
truncated after the second cumulant. This yields a gaussian distribution similiar to the Ornstcin-Uhlenbeck process. 
More recently, an explicit derivation of the Feyman path integral representation for the propagator of the radiative 
transfer equation has been given @. Here it was again evaluated by truncating the cumulant expansion after the 
second term. This was justified by declaring that photons are massless and non-interacting. However, the imposition 
of the speed constraint would not allow this gaussian approximation. 

In this paper, we describe the light propagation in stochastic media as the motion of a kind of Brownian particle 
on which the fluctuating forces act only perpendicular to the direction of its velocity. This is effective in strictly and 
dynamically preserving the speed of the particle. This process is shown to correspond to a diffusion in the angular co- 
ordinate in the velocity space for a white noise disorder. Exact expressions for the moments of the space variables are 
presented and the second cumulant approximation is shown to yield a gaussian expression similiar to the traditional 
Ornstein-Uhlenbeck theory of Brownian motion. An expression is derived for the probability distribution for large 
force strengths which preserves the light cone. The exact Fokker-Planck equation for the probability distribution is 
derived from the stochastic Langevin equations for a white noise process. Numerical solutions of this equation are 
presented . It is shown that the probability distribution in infinite media is strongly forward peaked for short times 
and randomizes only at times of about 8t* to lOi*. We have also solved numerically the equation for a semi- infinite 
geometry and obtained the persistence exponent of 0.435i ± 0.005 in 2-dimensions for this process. Solutions for 
a finite geometry are also given, showing that the effective diffusion coefficient as measured in a pulse transmission 
experiment through very thin slabs [L ~ I*) would be lowered. The effect of light amplification in the slab is examined 
briefly. 



II. THE MODIFIED ORNSTEIN UHLENBECK PROCESS 



Light scattering in a stochastic medium is treated as a probabilistic process where each scattering event only changes 
the direction of the photon. The wave nature and polarization effects are ignored and light is treated as a particle in 
a medium which exerts transverse fluctuating forces on the particle. It should be remarked here that while the actual 
disorder is maybe in space (quenched disorder), all current treatments including ours, are in terms of a Brownian 
motion (temporal disorder i.e., a stochastic process). This is a valid approximation for incoherent transport in the 
weak scattering limit (kl* >> 1 where k = A being the wavelength of light in the medium. The equation for 
the motion of a randomly accelerated particle with the special condition that the random forces always act only 
perpendicular to the velocity can be written as, 

?=7xf(t) (1) 
This we term as the modified Ornstein-Uhlenbeck process. We will consider two dimesions for simplicity, and write 

x = -f(t)y (2) 
y = f(t)x (3) 

where the force term f(t) is a random funtion of time. We will assume a delta-correlated force with gaussian 
distribution, i.e., 

(f(t)) = (4) 

(f(t)f(t'))=m-t f ) (5) 

and all higher moments of f(t) being zero. This makes our treatment most valid for a very dense collection of highly 
forward scattering weak anisotropic scatterers. This set of stochastic Langevin equations yield on integration a first 
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constant of integration x 2 + y 2 — c 2 , where c is the constant speed. So we can choose x = ccos9(t) and y = csin#(i) 
where 9(t) is some function of t. 9 is recognised to be the angular co-ordinate in the velocity space. Substituting 
these expressions back into equation(||,||), we obtain = f(t) or 

e(t)-e = f f(t)dt (6) 

Jo 

Hence 9(t) follows a Wiener process and we can write the probability for 9{t) as 
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This is the result for a diffusion in 9 the angular co-ordinate in the velocity space, and we recognise this modified 
O-U process to be a random walk on the circle of radius c in the velocity space. Constraining 9 to the range [0,27r], 
we get the marginal probability distribution for 9. 

n= — oo v ' v J 

The value of 9q can be conveniently chosen to be zero. 

Now we will derive the probability distribution function in the phase space. Consider the system of Stochastic 
Langevin equations 

x = u (9) 

y = v (10) 

u = -f(t)v (11) 

v = f(t)u (12) 

Let H(x,y,u,v) be the phase space density of points for the given system and U be the vector (x,y,u,v). Now, II 
satisfies the stochastic Liouville equation. 

^ + Vu • (tin) = o (13) 

where Vu = (J|, g^, -iz, -§^)- Substituting for U and averaging over all possible configurations of disorder, by the 
van Kampen Lemma |17| , the probability distribution P(x,y,u,v) = {H(x,y,u,v)) and satisfies, 

dP dP dP d d 

a + "15: + "W-^ </Wn> + "a; (/WnH0 (14) 

By the Novikov theorem jl8| for a white noise process f(t), 

</(W(t)]> = T -( ™ ) 

dP dP 

du dv 

Using the above, we obtain for P(x,y,u,v) the differential equation 

dP dP dP T / d d K ~ 



(15) 



at 1 u ^ + v ^^2 [ u d^- v d^) p (16) 



Now expressing u and v in terms of the angular co-ordinate 0, we finally get 



dP n dP . BP d 2 P 

- +C cos*- +C smtf- = 2r ^ (17) 
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This differential equation explicitly preserves the constancy of the speed of the photon. This Fokker-Planck equation 
is the same equation (in two dimensions) which was written down in Ref. [fL4f . It is rigorously proved therein that 
this has a path integral solution and the two approaches are equivalent. It appears that this equation has solutions 
in terms of Mathieu functions. However we have not been able to analytically solve the equation. 

The moments of the displacements can however be calculated analytically. The displacements can be written in 
terms of 9 as 



x — Xq 



y - yo 



c cos OdO 



sin 6d6 



Using these and a gaussian distribution for f(t), we get 

2c 
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(18) 
(19) 

(20) 
(21) 

(22) 

(23) 
(24) 



This reproduces the result of the traditional Ornstein-Uhlenbeck process in that the first moment saturates at a mean 
free path I* and the second moment increases linearly with time at long times (Tt/2 3> 1). For short times (rt/2 <C 1), 
the longitudinal spread (Ax 2 ) ~ t 2 and the lateral spread (Ay 2 ) ~ t 3 which are considerably slower than the diffusive 
linear behaviour. From these relations, we identify the mean free time t* to be p- and the transport mean free path 

2 

I* = ct*. The diffusion coefficient is identified as the coefficient of the linear term of the second moment i.e., 4r 



It is of interest to note that an analytic expression for moments of all orders for the displacements can be obtained. 
This expression is given in the Appendix. The marginal probability distribution function P(x, y, t; Xq, yo, 0) can be 
written in terms of a cumulant expansion (See the Appendix). Truncation of the cumulant series after the second 
term yields the result of Ref. jlij for the probability distribution. 



P{x,y,t;x o ,y o ,0) 



1 



Mr 



2ndet(M) 
2c 



exp ■ 



-(f - ro - a)i(r- ro - a)j 



(25) 



3= -(l-e- r */ 2 )(cos0 o ,sin# o ) 



r 



M tj = {(r- fo)i(r- fo)j ) - ( (r - f ) t )( (r- r )j ) 

The distribution is gaussian in this approximation and similiar to the distribution for the traditional O-U process 
fL2f . Thus it does not exactly preserve the light cone and would appear to constrain the speed only in an average 
sense. Higher cumulants would be required to describe this feature of fixed speed. 

An approximate solution which preserves the light cone can be obtained under the assumption that 9 is completely 
randomized in time t* , so that 6 has a uniform distribution over [0, 2tt]. This can be justified in the limit of large force 
strength (T), when the scattering events change the momentum by a large amount. Now the time can be discretized 
on this time-scale and the probability distribution can be written as 



P(x,t;x o ,0) 



1 

2^ 



duj e^-^i 



exp[ 



-iu>c'y* j cos ( 
j=i 



(26) 



where we have used that at t = 0, the angle 6 was uniformly distributed. To evaluate the average, we will use the 
fact that in this approximation, each 9j is independent of all others, giving 
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( exp[-ia;c^cos0,-i*] } = [J (uct*)] n (27) 

3=1 

where J is the ordinary Bessel's function of order zero. Using n = t/t* , we have 

1 f°° 

P(x, t;x ,0) = — du^'-^ [J^wd*)]*'* (28) 

Using the fact that the Fourier transform of J (ojct*) is zero for \x — xq\ > ct* and the fact that P(x, t, Xq, 0) is an n th 
convolution of Jo (wet*), it is seen that P{x, t, Xq, 0) is zero for \x — x Q \ > net* — ct. Thus the light cone is preserved. 
In fig- 1 , we plot the P(x, t, xo, 0) obtained by numerical evaluation for different times. It is seen that for t = t* , the 
probability is accumulated at the light front, and all the curves show a cut-off at \x — xq\ = ct. At long times, using 
the Laplace approximation, we have (for large n); 



[J${ujct*)] n ~ exp 



c 2 t* 2 uj 2 n 
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I \ 1/2 f u-.r>,) 2 



P(,,t;x ,0).(^ exp (29) 



Thus we recover the diffusion limit at long times. 



III. NUMERICAL SOLUTIONS AND RESULTS 

In this section , numerical solutions for the differential equation ( |i"7j ) are presented. The particle is released in the 
x — y plane at the origin (generally) along an initial direction 9q. Here 9 is the angle made by the velocity vector 
with the x-axis. Let us first further simplify by assuming invarience with respect to y i.e., we have a line source 
along the y axis. Then the derivative with respect to y drops out and and we have a partial differential equation 
in three variables. This is essentially a parabolic equation with an advective term. To numerically propagate the 
probability distribution in time, we use an alternating direction implicit -explicit method ]l9[ | for x and 0. A local von 
Neumann stability analysis [jl9f shows that this differencing scheme is unconditionally stable. The initial condition is 
a (^-function at x — 0, 9 — which is approximated by a sharp gaussian for numerical purposes. For infinite media, 
the boundary condition P(x, t) = for |x| > ct is used. For a semi-infinite medium — oo < x < L with an absorbing 
boundary at x = L, the appropriate boundary condition is given by P(L,8,t;xo,9a,0) = for — n < 9 < — 7r/2 and 
7r/2 < 9 < 7r, corresponding to no flux entering the medium from free space. Also, we can write the Fokker-Plank 
equation in the form of a continuity equation. 

£>P 

7* +v -' = ° (30) 



v = ex d-x +ee de 

dP 

j = e x cos9P + e 9 T — 
o9 



Since T = outside the medium, we can conclude that the current density J* in the real (x) space is conserved across 
the boundary in the forward direction (—tt/2 < 9 < ir/2) while the current density in the velocity (9) space is not 
conserved. This explains why the output flux at the boundary is proportional to the value of the probability distri- 
bution function at the boundary itself ( rather than the space derivative of the probability distribution (§7) given 
by Fick's law ) as observed in experiments f20[| . For a finite slab we use a similiar boundary condition at the other 
boundary. 

In Fig-2, we show the probability distributions in an infinite medium with the initial condition, P(x,9,t = 0) = 
5(x)5(9). It is clearly seen that the probability distribution for times upto 5t* is peaked in the forward direction 9^0 
for x > 0, with a tail in the backward direction (9 ~ ±7r) at x < 0. There is also a clear cut off at |x| = ct, which 
is prominently noticeable for positive x. The small amount of tailing arises from the finite width of the gaussian by 
which the (5-function was approximated. One can also note that the probability distribution becomes almost flat along 
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the 0-axis only at times of about eight times the mean free time (8t*). In Fig-3, the first and second moments of the x 
co-ordinate are shown. The solid lines show the analytical results of equation (^(] and |2^) and the symbols(o) represent 
the results of the numerical solutions. Excellent agreement is found between them. In Fig-4, we show the marginal 
probability distribution for x i.e., P(x, t; xq, 9q, 0) = J_ d6P{x, 9, t; xq, 9q, 0). At short times (t ~ 3i*), there is a clear 
ballistic peak, separate from the more randomized tail. The probability distribution for these times is also clearly 
forward-peaked. One can also note that the probability distribution randomizes and becomes almost gaussian, centred 
at x ~ /* only at times t > 8t*. As noted above, this is also the time by which the angular coordinate 9 randomizes. 
This is when the diffusion approximation becomes valid. This can be understood by noting that, by equation^, 
the time required for Pt{9) to attain an angular width of 2ir is T where T is given by (A9 2 ) — (2tt) 2 ~ 2TT. This 
yields (using Y/2 — t*) a value of T = ir 2 t* ~ lOt* for the randomization time. Thus we now have a clear picture of 
the reason for this long known experimental fact ||. This forward peaked behaviour at short times also illustrates 
the deficiency of the second cumulant approximation where the probability distribution is a Gaussian and symmetric 
about the first moment. Higher cumulants are clearly required to describe these asymmetric features. 

The probability distribution functions for a semi-infinite medium are shown in Fig-5. Here the particle is released 
at the origin inside the random medium and the initial direction is towards the boundary (in this case at x = Al*) For 
times lesser than At* , there is no difference in the probability distribution from the case of the infinite medium. This 
is because the wave front has not propagated upto the boundary and the effect of the boundary is not felt. This is 
to be contrasted with the diffusion approximation where the effect of the boundary is felt everywhere simultaneously 
and causality is violated. At long times the probability distributions attain a typical shape with a long tail at negative 
x within the medium and a sharp cut-off at the boundary. In Fig-6, we show the marginal probability distribution 
for x i.e., P(x,t;xo,9 ,0) = J^d9P(x,9,t;xQ,9 ,0). The value of P(x, t; Xq, , 0) is finite at the boundary and 
zero outside. As seen in Fig-6b, if the points near the boundary are linearly extrapolated outside the boundary, 
they all roughly cross the x-axis at about 0.72* which is the value of the extrapolation length used in the diffusion 
approximation HJ. In Fig-7, the surviving probability inside the medium P s = J dx J d9P(x, 9, t; xo, do, 0) is plotted 
with time. For long times, this quantity should scale as f - * where $ is the persistence exponent for this process [ p2| . 
We have performed these calculations for several source-boundary distances and obtained a value of 0.435 ± 0.005 as 
the persistence exponent for this process in two dimensions. 

Finally we present solutions for a finite slab with absorbing boundaries at x — ±L. The particle is released from 
the origin at t = along the positive x direction. Fig-8a shows the first and second moments of the probability with 
time in a thin slab of thickness 21*. The first and second moments initially increase as in an unbounded medium 
until the photon-front hits the boundary and dip before increasing again and saturating at an almost constant value. 
The dips occur because just after the ballistic and near ballistic components exit the slab, only the photons which 
are effectively moving in the opposite directions are left behind. In fact, the first moment is seen to become negative, 
implying that the net transport is in the backward direction for some time. The dip in the second moment implies that 
the photon cloud is effectively expanding at a slower rate. This would cause a lowered 'effective diffusion coefficient' 
to be measured in a pulse transmission measurement. This reinforces the conclusions reached in Ref. Jl3j based on 
Monte-carlo simulations and explains the experimental results of Ref. || on a more rigorous footing. Fig. 8b shows the 
survival probability for the case of a finite slab. This decays considerably faster than the in case of the semi-infinite 
slab, though at early times (t ~ t*) the decay rates are comparable. The initial rates of decay are comparable because 
of the forward peaked nature of the probability distribution at early times, when the effect of the boundary at the 
back is hardly felt. This is to be compared with the mirror-image solution in the diffusion approximation, where equal 
weightage is given to both boundaries at all times. 

Finally we turn to the case of an amplifying stochastic medium. The effect of medium gain can be incorporated 
straight forwardly by noting that in our treatment, the time of exit from the slab directly translates into path-length 
traversed within the medium because speed is kept absolutely fixed. In the presence of amplification in the medium 
therefore, the net gain is directly proportional to the time. Thus the output flux at the boundary in a given direction 
is simply P(L, 9, t) cos 9 exp(at), where a is the gain coefficient in the medium. It is thus simple to obtain a picture of 
amplified emission from such a medium. In Fig-9, we show the total light emitted by slab with boundaries at x — ±21* 
for several amplification factors. The photon is released from the origin in the positive x - direction. For large times, 
the output increases exponentially because of the presence of a exponential gain in the medium with no saturation. It 
is seen that the ballistic part is only slightly amplified while the output in the tail regions are increased considerably. 
To obtain a more realistic picture of lasing in random media p3| , p4 however, one would have to consider the lasing 
level population depletion and saturation effects. 
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APPENDIX : EXPRESSION FOR THE MOMENTS (X N ) 
The n th order moment is given by 

((x - x ) n ) = c" f [ ■■■ f dtndtn^f-dhicoseitx) cos 9{t 2 ) ■■ -cos 6(t n )) 
Jo Jo Jo 

Writing 6(ti) as di, the quantity within the angular brackets can be expressed as follows, 

(cose*! cos0 2 ---cos0 n > = 2~ n < (e 101 + e^ 1 )^ 2 + e"^ 2 ) • • ■ {e lB - + e~ lB -) ) 



= 2-" Y. < exp[*5>i0j] ) 

3=1 



cr;=±l 



This can be expressed as a path integral using a gaussian distribution for f(t) 



exp[z °M ) = / W)} ex P | - / 



2r 



3=1 



dt' 



n I n 



exp 



where to = and we assumed a time ordering of t\ < t 2 < • • • < t n . Thus, 



t ft„ 
o Jo 



n I n 



{{x-x ) n )=c n (nl)2- n / dt n / *„-!••• / dii ffi P -^E E"i (tfc-rjk-x) 



T1 ,CT 2 ■■ -CTt, 

<Ti=±l 



fc— 1 Vj — k 

A similiar expression can be obtained for the ((y ~ yo) n ) by noting that sin# = cos(7r/2 — 9). 
Now we can obtain the joint probability distribution of x and y as 

P(x,y,t;xo,yo,0) = { 5(x — x — c cos 9 (t')dt') 8{y — yo — c / sm9(t')dt') ) 



exp 



Expressing the ^-functions in terms of the Fourier transforms, 

2 



P(x,y,t;x ,y ,Q) 



duj x dujy e *K(z-zo)+^(y-w)]( 



(31) 



(32) 



(33) 



(34) 



-ic / [oj x cos9(t') +U! y sm9(t')]dt' 
Jo 



(35) 

) 

(36) 



This statistical average can be evaluated by a cumulant expansion |15| and since we have an expression for moments 
of all orders, we can in principle evaluate the cumulant expansion to any desired order. 
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FIGURE CAPTIONS 



Figure-1 : The marginal probability distributions P(x, t; xq, 0) predicted by the approximate solution given by equa- 
tion(28) at different times indicated in the figure. There is a clear cut-off at the light front and initially the probability 
accumulates at the light front (for t = t*). 

Figure-2 : The Probability distributions in the phase space of a particle in an infinite medium at different times ob- 
tained by numerically propagating equation(17). The particle is released at x — along the positive x direction^ = 0) 
at t = 0. The probability distribution is clearly forward peaked and and becomes almost flat along the #-axis only at 
times of about 8t* . 

Figure-3 : The first and second moments of the displacement for the probability distribution of a particle in an infinite 
medium. The solid lines show the anlytical result of cquation(20) and equation(22) while the symbol ( o ) show the 
result obtained from the numerical solutions. 

Figure-4 : The marginal probability distribution P(x, t; Xo, #o> 0) = P(x, 9, t; xo&o, fydO at different times. The 
marginal probability distribution becomes almost a Gaussian at times of 8t* . 

Figure-5 : The Probability distributions in the phase space of a particle in a semi-infinite medium at different times. 
The particle is released at x = along the positive x direction(# = 0) at t = 0. The absorbing boundary is located 
at Al*. The probability distribution is zero in the range — 7r < 9 < —ir/2 and ir/2 < 8 < ir at the boundary implying 
there is no incoming flux into the medium. 

Figure-6 : The marginal probability distributions in a semi- infinite medium with an absorbing boundary at x = Al*. 
The plot on the right shows an expanded view of the distributions near the boundary. The solid straight lines are the 
linear extrapolations of the behaviour near the boundary. All of them are seen to cross the x-axis roughly at 0.7Z* 
outside the boundary. 

Figurc-7 : The surviving probability of the particle inside the semi-infinite medium for an absorbing boundary at 
4P(o) and 2P(*). The persistence exponent •& is obtained from the long time behaviour of the survival probability. 
The lines (• • •) and (— • — ) show the linear fits and give a persistence exponent of 0.4309 and 0.4364 respectively. 

Figure-8 : The first moments of the displacement of a particle in a finite slab of thickness 21* (left plot). The right 
plot shows the survival probability in a semi-infinite medium and a finite slab. The distance between the point where 
the particle is released and the boundary is same in both case. (2/*). 

Figure-9 : Total light emitted (from both sides) by a disordered slab with amplification for different values of the gain 
coefficient A in the medium. The output increases exponentially at long times. 



9 




1 2 3 2 4 

X (in units of I*) X (in units of I*) 




0.2 




X (in units of I*) 



t = 5t* 

t = 3t* 




0.04 



t = 7t* 



0.01 



t = 32t* 




1 2 3 4 5 



log(Time) (Time in units of t*) 



